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Abstract 

We reanalyze the contributions to neutrinoless double beta (fiPov) 
decay from supersymmetry with explicit breaking of R-parity. Although 
we keep both bilinear and trilinear terms, our emphasis is put on bi- 
linear R-parity breaking terms, because these mimic more closely the 
models where the breaking of R-parity is spontaneous. Comparing the 
relevant Feynman diagrams we conclude that the usual mass mechanism 
of double beta decay is the dominant one. From the non-observation of 
PPov decay we set limits on the bilinear R-parity breaking parameters 
of typically a (few) 100 keV. Despite such stringent bounds, we stress 
that the magnitude of R-parity violating phenomena that can be ex- 
pected at accelerator experiments can be very large, since they involve 
mainly the third generation, while PPq u decay constrains only the first 
generation couplings. 
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1 Introduction 



The most well-studied realization of supersymmetry has a conserved R-parity. 
This is a multiplicative quantum number which can be defined as Rp = 
(— l) 3B+i+2 ' s , where B and L are the baryon and lepton number and S the spin 
of the corresponding particle. This defining property of the Minimal Super- 
symmetric Standard Model (MSSM) ffl] is theoretically ad hoc since the origin 
of R-parity conservation is unknown. Since the early days it was recognized 
that the violation of R-parity could arise as residual effects in larger theories |2| 
or spontaneously due to a sneutrino vacuum expectation value ( VEV) H . We 
find the models with pure spontaneous breaking more well-motivated than 
explicit breaking models, in the sense that they put R-parity breaking on a 
dynamical footing. They have only lepton number violating ) terms and 
are therefore in agreement with proton stability. The simplest models of this 
type (with the SU{2) ® U{1) gauge structure) are characterized by the exis- 
tence of a physical Goldstone boson, generically called majoron. In the early 
versions the majoron was an isodoublet (the imaginary part of the sneutrino) 
[[J; these are now ruled out by LEP measurements of the invisible Z width [Q. 
However, in the modern variants of supersymmetric models with spontaneous 
R-parity breaking the majoron is a singlet (such as a right-handed sneutrino) 
and therefore can not be produced in the decays of the Z. Such {$, p ) models are 
both theoretically attractive and phenomenologically viable [|5], However, 
quite generally study of broken R-parity phenomenology has received a lot of 
attention recently 0, || |9|, [T0|] . 

Especially interesting in view of the mounting evidence for non-zero neu- 
trino masses |ll|, 112, |H| is the fact that Jji p models allow for finite - and in 



principle calculable - neutrino masses, and moreover could even explain their 
observed smallness without invoking some high-energy scale physics |14 |, like 
one does for example in the see-saw mechanism 

Neutrinoless double beta (P/3q u ) decay is a AL = 2 process and therefore 
one naturally expects it to occur in models with lepton number violation in 
the Lagrangian. Even though there is a variety of mechanisms inducing (313q v 
decay in gauge theories, one can show that whatever the leading mechanism 
is at least one of the neutrinos will be a Majorana particle |]I6 ], as illustrated 
in the black-box diagram of Fig. 1. This well-known argument establishes 
a deep connection between Majorana neutrino masses and /3/3o u decay: in 
gauge theories one can not occur without the other being present. The same 
remains true in supersymmetric theories |T7|], where moreover one can show 
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that also the supersymmetric partner of the neutrino must have a B — L 
violating Majorana-like mass term, if a Majorana mass of the neutrino exists 
[[j] ]. Turning the argument around, one expects that the observed absence of 
PPou decay allows to derive stringent limits on Ifl p parameters. This has been 



shown for models with explicit trilinear R-parity breaking in |i9fl . In the 
current work we will concentrate on the bilinear lji p terms. 

/3/3q u decay within supersymmetric models with bilinear R-parity breaking 
has been considered previously in f20| . In our derivation the final expression 



of the decay rate (and the corresponding nuclear matrix elements), however, 
differ from the one given in | 20| , details are given in section 3. 

In an exactly supersymmetric and S , [/(2)i-invariant world R-parity Vio- 
lation (RPV) can be moved around the Lagrangian by field redefinitions [Q. 
This observation lead to the claim, that the bilinear RPV terms are unphysical 
and should be discarded. However it has been shown that in realistic softly 
broken supersymmetric models it is in general impossible to eliminate the bi- 
linear terms [plfl . Moreover, as already stated, it is possible that R-parity is a 
symmetry (maybe approximate) at the Lagrangian level, broken only by the 
ground state f5|. Such scenarios lead only to fj being therefore consistent 
with baryon number conservation. Models with bilinear terms can be imag- 
ined as a truncated version of the full dynamical models |2l| and as such offer 
a systematic way to investigate RPV at the phenomenological level. 

This paper is organized as follows. In the next section we introduce our 
conventions and notations. Section 3 derives the effective Lagrangian of (5(5^ 
decay in ^? P SUSY with bilinear terms, while section 4 summarizes the necessary 
nuclear physics input. Section 5 then analyzes the constraints P/3 0u decay can 
put on bilinear lfi p parameters. 



2 Super symmetry with explicit R-parity break- 
ing terms 

In this section, we set up our definitions and notation. Our conventions are 
usually such that in the limit of vanishing I/l p parameters the MSSM notation 
of H is recovered. 

In the presence of bilinear R-parity and lepton number violation there is no 
distinction between the lepton doublet and the Higgs doublet superfield giving 
mass to the down-type quarks. This fact can be accounted for by defining a 
superfield $ as 
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$ = {H 1 ,L 1 ,L 2 ,L 3 ). 



(1) 



For the MSSM field content the most general gauge invariant form of the 
renormalizable superpotential can then be written as 



W = e ab [^Vf&jEg + \T® a iQ)D c k + K k Q«H b 2 UC + ^^H b 2 \ . (2) 

Here, Q and D c , U c are the quark doublet and singlets superfields, respec- 
tively, E c is the lepton singlet superfield and H 2 the Higgs superfields with 
Y(H 2 ) = 1 responsible for the up- type quark masses, with h? k being the cor- 
responding Yukawa couplings. The indices j,k = 1,2,3 denote generations, 
whereas I, J = 0, 1, 2, 3. The indices a, b are SU(2) indices. In the basis (P 
one can separate W into an R-parity conserving and an R-parity violating {IfL p 
) part 

W = W Rp + Wfr, (3) 

where the Rp part is simply given by taking I or J to be zero and i,j respec- 
tively. Then X°J k = |/i{ fc and X^ k = . The Iflp part contains the remaining 
components, 

W % = \fUL 3 E c k + XfUQM + e t L t H 2 (4) 



where //j — > e^, to make contact with the notation of If W$ p would be 
the only source of RPV one could easily rotate the field $ into a basis $' with 
fi 1 = (fi', 0, 0, 0) effectively eliminating the bilinear terms. However, another 
source of RPV is found in the soft supersymmetry breaking part of the scalar 
potential. It contains the terms: 



V£ ft = Ai Jk <f>i$jEk +A I d jk $ I Q j D% + B I $ I H 2 

+ {m 2 IJ + f i I (ij)$ I &j + --- + h.c., (5) 



where the dots represent terms not interesting for the discussion here. Rotating 
the superpotential as discussed above, it is easy to see that as long as the 
B 1 are not exactly parallel to the /i / the effects of the bilinear terms of the 
superpotential in the rotated basis will reappear in the soft SUSY breaking 
terms. 

Let us briefly discuss the appearance of bilinear terms B I ^jH 2 . The pres- 
ence of these terms imply that the tadpole equations for the sneutrinos are 
non-trivial, i.e. the VEVs of the sneutrino fields in general are non-zero |E2|. 
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As a result the MSSM vertices Zvv and Wee create gaugino-lepton mixing 
mass terms Zv < v > and We < v >. This mixing then leads in the physi- 
cal mass eigenstate basis to lepton number violating (LNV) interactions. The 
LNV vertices relevant to (3(3q v decay will be discussed in the next section. The 
mass eigenstate fields are defined from the weak eigenstates by means of an 
orthogonal transformation, 

= Hi^'oj, (6) 
*(±)i = Ag*' (±)j -, (7) 
with the weak basis states in two component notation defined as 

% T = ^l^l^l,^',^^^), (8) 
*£) = (Vl,^ 2 ,Vi,-*A_,V>U ( 9 ) 

^ = (^,^,^3,-^,^), (10) 

where ^>\. (tpl ) are the neutrino (charged lepton) fields, and — i\ 3 ,\^ 

1 2 

are the unmixed photino and gaugino states, respectively, and tpjj are the 
neutral components of the Higgs doublets, and the Higgs and lepton fields are 
written in the original basis $. 

The mixing matrices S and A diagonalize the neutrino-neutralino and 
charged lepton-chargino mass matrices. Details are given in the appendix. We 
note in passing that the neutral mass matrix has such a texture that only 
one neutrino gains a tree-level mass, even if all of the e« are non-zero. The 
degeneracy of the remaining two neutrinos, massless at tree-level, is however 



lifted at the 1-loop level |23] 



3 Effective Lagrangian of (3(3q u decay 

Here we derive the effective Lagrangian describing (3(3q v decay in the RPV 
model. Formally one can define the amplitude for /3p 0u decay as, 

R Pfhv = ((A,Z + 2),2e-\T[exp(t J d 4 xC eff (x))]\(Z, A)}. (11) 

For the derivation of £ e // h is convenient to decompose the different I/ip 
contributions as shown in figure 2. As indicated in the figure, in the Jjt p MSSM 
the effective lepton number violating j3j3 Qu decay amplitude arises either by the 
AL = 2 (massive) Majorana neutrino propagator (fig. 2. a), or through short- 
distance four- and six-fermion lepton number violating vertices (fig. 2.b,c). 
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The latter appear in this model in two ways, either directly from the trilinear 
parameters of the superpotential eq. ([|), or indirectly from the bilinear terms, 
due to the neutrino-neutralino and char gino- charged lepton mixing. 

From the superpotential one can read off the following AL = 1 trilinear 
terms which are relevant for /3/3q v decay, 



C 



(Trilinear) 



(Af 



h>iP R d k )d* L - (ujP R d k )e* L - (eiP R d k )u* j L 
iP R e?)d k , R } + \f(v?P L d 3 )dl R + h.c. (12) 



As discussed above, the bilinear terms induce mixing between the SM leptons 
and the MSSM charginos and neutralinos in the mass-eigenstate basis, where 
one finds the following AL = 1 terms Ejfl, 



C 



(Bilinear) 



9_ 

V2 



^- Kn W-(e^P LX n ) + 



V2g {^(u k P R u c )u L + (3i{V k P R d)d R + fl e ki (p k P R e c )e u 
+ (3 c {uP R e c )d L ) + h.c. (13) 



The first term in eq. (|13[) is generated from the Standard Model PF-boson- 
fermion-fermion and R-parity-conserving interactions of the type W — x° — 
present in the MSSM. The R-parity-violating terms arise from the x ± -fermion- 
sfermion and x°-fermion-sfermion interactions. 

The coefficients in eq. (|13|) are defined as follows, 



K, 



k=l 



V2A7 d E 



— ^* 

14^n+3,5 



A - ^* 
^15^1+3,6 



tan^SfcA 



6 



(tan 9w^k,4 + 3Sfc j5 ) 



Indices i, k above run from 1 to 3, whereas the index n = 1,2, 3, 4. 



(14) 

(15) 
(16) 
(17) 

(18) 
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If the ftp parameters are smaller than the other typical supersymmetric 
parameters, one can find an approximate analytic diagonalization for the ro- 
tation matrices H and A using the method introduced in ref. P5|. This was 



done in ref. j26| and the results are presented in the appendix. In the analysis 
presented later on we will use this convenient approximative method and jus- 
tify its use a posteriori. In leading order in the expansion parameters £ and 
£l the coefficients defined above are given by, 

«» « £ K m C lm - V2^ L ) n N* n2 - (6O12K3 (19) 

m=l 

/^^tan%£(K*); fe £*i (20) 

P'k « \ E(V;) jk (tane w ^ + 3C J2 ) (21) 
3=1 

Pu « ^l)xi{V:)u - \ ETOifcCtan^j + (22) 

(23) 

In the derivation we have always neglected terms which are of higher power in 
the neutrino mass and terms which are suppressed by factors of m Ut d, e /M S usY 
due to left-right sfermion mixing and/or higgsino couplings. Here, m u ^ e are 
the masses of the up, down quarks and the electron mass and Msusy is the 
typical mass scale of SUSY particles, assumed to be of order of 0(100 GeV). 

In addition to the lepton number violating Lagrangians eq. (|12"D and eq. 
( |i"3"D for the calculation of C e ff for the (3(3q v decay the ordinary Standard Model 
(SM) charged-current interaction as well as the MSSM gluino-quark-squark 
interaction 

x(o) 

C~ g = -V2g 3 ^(q a P R ~g^f L - q a P L g^f R ) + h.c. (24) 

are needed. Here X^l are Gell-Mann matrices (a = 1,...,8), and the super- 
scripts a and j3 are color indices. 
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The leading order diagrams are shown in fig. 3. The first diagram, fig. 
3. a is nothing but the well-known mass mechanism, supplemented by a con- 
tribution from neutralino exchange []. Since for finite values of the bilinear 
parameters one neutrino acquires a tree-level mass this contribution is always 
present in the bilinear l/l p model. The diagram fig. 3.b represents a short-range 
contribution to (3(3q v decay induced by the non-zero chargino-electron mixing. 
The remaining two diagrams are of long-range type, and involve the standard 
Dirac-type (AL = 0) neutrino propagator. Note, that using the above approx- 
imations these diagrams are always proportional to the product of a bilinear 
times a trilinear parameter. However, these diagrams are also inherent to the 
model, even in the basis where the trilinear parameters are set to zero. Their 
contributions are, however, suppressed by a factor of mf/mj from left-right 
sfermion mixing, which should naturally be very small for the first generation. 
The diagrams in fig. 3 do not include those proportional to (A^ fc ) 2 . These 
have been previously analyzed in 



Integrating out the heavy fields and carrying out a Fierz transformation, 
one obtains finally the effective Lagrangian 

C eff (x) = V2G F {V v ) lk {e lli P L v k )J" 

+ ^{m + V x )J"J,{eP R e c )- (25) 
- V2G F (u k P R e c ) J + 77$ {vk^P R e c ) 
The first term in the above equation is the usual SM charged- current interac- 



tion. Color-singlet hadronic currents in eq. (|25|) are given by 



J = u a P R d a , (26) 

r = u a lt ,P L d Q} (27) 

= u a a ^P R d Q . (28) 

In eq. (^5|) we have introduced effective lepton number violating parameters 77 
defined by, 

1 Somewhat arbitrarily we call the 4 heavy neutral states neutralinos, while the light states 
are referred to as neutrinos. This is done, since in the limit of vanishing parameters the 
MSSM neutrino and neutralino states are recovered. 
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(30) 

*7a =7^-1 2 — —Oil H 2~Oifef (31) 

/(tip >- mi mi mi mi > 



Cifl 9 9 "'6-L 9 1 9 

zGf l mj m~ m~ wA 

e L,i U L dji 

(k) _ g^d 1 r Pk x Pi x , P 



{— ^-Sii H 2~kk\ (32) 

l. mi mi mi > 



q 8Gf ^rnj,/ lL m\ ~ tx ' m~ 

U L d R d L 



Here, Gp is the Fermi constant and m p is the proton mass. 

At this point, we would like to discuss some differences between our cal- 
culation and an earlier work [20Q. In ref. [^0| it was found that the short- 



ranged part due to gluino exchange is proportional to pseudo-scalar and ten- 
sor hadronic currents. (This would imply an important enhancement of the 
gluino diagram due to the pion exchange mechanism first discussed for the 
case of trilinear R-parity breaking in fll9| .) However, we find that the gluino- 
mediated diagram is proportional to the same hadronic structure as appears in 
the heavy Majorana neutrino exchange, characteristic of left-right symmetric 
models |37L |23j. Therefore, we conclude that the gluino contribution is negli- 
gibly small. Also in the long-range part of the effective Lagrangian, eq. (f2~5|), 



the last term in eq. ( |3"TD is missing in the corresponding equation in ref. [EU 
and the long-range tensor contribution proportional to r)^ has been neglected 



compared to the long-range scalar-pseudoscalar part in [20 . 

The effective Lagrangian given above is defined in terms of quark currents. 
For /3/3 0u decay we have to reformulate £ e // in terms of nucleon currents and 
apply the proper nuclear wave functions to describe nucleon states in the nu- 
cleus. This is done in the following section. 



4 Half- life of /3/3q u decay and definitions of nu- 
clear matrix elements 

Having specified the effective Lagrangian eq. ( p5| ) the derivation of the half- 
life of double beta decay is rather straightforward. Following the standard 
procedure described in [Tj| |30[ , after lengthy manipulations one finds, 
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V 1/2 



-1 



G 01 (M v y\v 



2 1 AL=2|2 



(33) 



Here, Goi is the leptonic phase space integral, numerical values can be found 
in ||27J, Ai u is the nuclear matrix element governing the well-known mass 
mechanism of double beta decay and 



V 



AL=2 



(m„) . TfLp i r ry i 



where T^- = M.j/M. v with j = A, Ag, N. 

Individual nuclear matrix elements are defined by, 



M u = 
M N = 



(M Gr , m -(^) 2 M F , m ), 
9 A 

(M GT . N -(^) 2 M FtN ), 

9A 

= a 1 (-M GT '+M T ,), 



M X ,q 

The coefficients a, are given by p9| 



( — aiM GT i + a 3 M F ' + a 2 M T i). 
3 



,(3) 



a i 



4Tj 3) g v (l - 2m p (g w /g v )) 
g\(Rm e ) 



(Rm e )g A ' 



a 2 



a 3 



4(2T 2 (3) -T± 3) )g v 



(34) 



(35) 
(36) 

(37) 
(38) 

(39) 
(40) 
(41) 



g\{Rm e ) 

with gA — 1.26, gy ~ 1 and (gw/dv) — — 3.7/ (2m p) from CVC, while numer- 
ical values for can be found in ref. ||31|| . m e and R are the electron mass 
and the nuclear radius. The matrix elements in the closure approximation are 
(summation over nucleon indices are suppressed) given by |32|, |29[ 



M GT)k = (0ph k (a a a b )T+r+\0t) 
M F;k = (0+|^r a + r+|0+) 



(42) 
(43) 
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M T > = (tf\h T ,{{a a ? ab ){<7 b ? a b) - ~(^&)K + t 6 + |0+). (44) 
with k = m, R, N and and h T > are neutrino potentials defined as 

h m = -R dqq — , 45 

71 Jo uj(oj + E) 

. 2i? 2 r» 4 j (gr afe )/ 2 (g 2 ) 

^ = / ; — nPV - ' 46 

TTmpJo uj(uj + E) 

V = -— / dqq 2 ^-^Lr{q 2 j G {qr ab ) - 3— ji(gr a6 )}. (47) 
limp Jo u(u + E) r ab 

= ±r ab / ^ e tqr °" 48 

m p m e J (27i)- i [mr A + q^) 4 



with xa = fTtA^ab, = y q 2 + m 2 , E the mean excitation of the intermediate 
nuclear states and r ab the nucleon-nucleon distance. f 2 (q 2 ) is the nucleon 
form factor, usually taken to have dipole form, jo and j± are spherical Bessel 



functions, and tua is typically ttla = 0.85 GeV f32fl . 



5 /3/3q u Decay Constraints on fi p Parameters 

Up to now our calculation has been independent of the use of any nuclear 
structure model. However, for a numerical analysis we have to assume some 
definite numbers for the relevant nuclear matrix elements. For defmiteness we 



will use the pn-QRPA values of the nuclear matrix elements QTSJ, [28], p2| . 
Note that the uncertainty of limits on the fundamental ft p parameters depend 
only on the square roots of the uncertainties of the nuclear structure matrix 
elements. 

In the following analysis we will use the currently most stringent limit on 
double beta decay 

T 1/2 { 76 Ge) > 1.1 x 10 25 ys (49) 



found by the Heidelberg- Moscow experiment [[33j . 

The experimental constraint eq. ( [49"D together with the theoretical decay 
rate eq. ([£J) defines an excluded area in a complex multi-dimensional pa- 
rameter space consisting of the ft p parameters as well as a number of MSSM 
parameters. 
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We have analyzed the contributions from the different Feynman diagrams, 
see fig. 3, and found that for typical choices of MSSM parameters the neutrino 
mass mechanism is always the dominant contribution. To understand this 
finding, let us first explicitly calculate the effective neutrino mass in the bilinear 
-^pHiodel, which, using the usual GUT assumption | M x = | tan 6^M 2 , is given 
by, 

where A x = E\V\ — LO\H (see appendix) and the prime indicates summation only 
over light fermion states (not neutralinos) , and u)\ := (vi). 

Let us first compare the short-range contributions with the long-range 
(light) neutrino exchange. All of them are of the same order in RPV pa- 
rameters, i.e. proportional to A^. By naively comparing the corresponding 
diagrams one might think that they are of equal importance. Closer inspec- 
tion, however, reveals that this is not the case, and that the short-range con- 
tributions are suppressed compared to the long-range ones. This suppression 
is due to nuclear physics effects. For heavy particle exchange the two quarks 
interacting have to come very close together, where the strong repulsive part 
of the nucleon-nucleon interaction is important. For point-like nucleons this 
contribution would be zero, so that the corresponding nuclear matrix elements 
become non-zero only due to the finite nucleon size. A naive estimate of the 
size of the short-range contributions relative to the neutrino mass contribu- 
tion can be given by comparing the typical momentum scale of the nucleons, 
order p F , with the typical mass scale of SUSY particles m S usY and thus is 
expected to be Pf/^susy ~ 1CT 3 . Numerically we have found that indeed 
the constraints on RPV parameters derived from the gluino diagram alone are 
weaker than the constraints derived from considering light neutrino exchange 
by about the above factor. 

In contrast, for the long-range diagrams, fig 3.c and 3.d, an enhancement 
compared to the neutrino mass diagram would be naively expected. The 
explicit calculation, however, shows that these diagrams are proportional to 
A^ 11 Ai (in a basis with non-zero trilinear terms) or proportional to rrif/mjAi ~ 
(9(10 < - _6 ^)A 1 for first generation fermions (in a basis where A^ 11 is assumed to 
be zero) . Thus despite a huge enhancement factor of these diagrams, from the 

2 Although we will use this assumption in the following, our results do not sensitively 
depend on this choice. 
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propagation of the neutrino a la Dirac, no useful limits on the fundamental 
parameters e\ and uo\ can be derived from these contributions. 

Fig. 4 then shows the limits on the parameters e\ (fig. 4. a) and u\ (fig. 
4.b), assuming tan/3 = 1 as a function of the MSSM parameters (M 2 ,/i). 
Limits typically of the order of a few hundred keV are found, choosing M 2 and 
fi of the order of up to a few hundred GeV. Note, however, that the limit on 
ei depends rather strongly on the choice of tan /3, and detoriates proportional 
to tan /3 for large values of this parameter. 

The robustness of the double beta decay bound is then demonstrated in 
fig. 5, where calculated half-lives for the (3f3o v decay of 7e Ge as a function of 
(b>i) are shown. For this figure we have varied the MSSM parameters freely 
and randomly in a range motivated by naturalness arguments, i.e. M 2 and /i 
between 100 GeV and 1 TeV and tan/3 = 1 — 50. A clear absolute upper limit 
of the order of 1.5 MeV on (pi) is found from this analysis. 

The limits shown in fig. 4, 5 assume that only one bilinear ]ji p parameter 
is non-zero, which is unnatural from the point of view of the minimization of 
the scalar potential ||2TJ. A combined limit on (ei,Ui) is therefore shown in 
fig. 6., for fM = 100 GeV, various values of M 2 and tan/5 = 1 (fig. 6. a) and 
tan/3 = 10 (fig. 6.b). Obviously, the limits on e\ depend rather strongly on 
tan/3, whereas the limits on uo\ do not. Note the double logarithmic scale. 

In summary we find that j3f3o v decay provides very stringent limits on RPV 
parameters of the first generation. For any reasonable choice of MSSM param- 
eters in the bilinear ^ p model the mass mechanism is the dominant contribution 
to /3j3 0u decay. 

6 Discussion and Conclusion 

We have calculated the contributions of (bilinear) R-parity breaking super- 
symmetry to neutrinoless double beta decay. This model mimics closely the 
case of spontaneous breaking of R-parity. We find that [3(3q u decay constrains 
only a subset of the possible bilinear parameters, namely e\ and the sneutrino 
VEV of the first generation u\. This is a general property of the theory and 
does not mean any fine-tuning of parameters. For the first generation lfi p 
parameters, on the other hand, Pf3 0l/ decay provides very stringent limits, typ- 
ically of the order of a few hundred keV up to a few MeV. With these limits 
at hand it seems rather hopeless to search for RPV in the first generation at 
colliders. From this point of view a possible interpretation of HERA events in 
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terms of ljl v interactions would be rather unlikely. However, it is important 
to stress that, even though we find rather stringent constraints on the mag- 
nitude of first generation ljl v parameters, these do not limit in any way the 
attainable magnitudes of R-parity breaking signatures expected at colliders, 
since the latter involve mainly the third generation, i.e. they involve r's or 
v T 's. This fact follows very naturally, for example, in scenarios with radiative 
symmetry breaking, due to the larger Yukawa couplings that characterize the 
third generation. 

In our present study we find that the main origin (3(3q v decay in this model 
is the mass mechanism and conclude that other contributions are practically 
irrelevant. 
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7 Appendix A: Mass matrices 
7.1 Neutral fermion mass matrix 

The presence of the lepton number violating bilinear terms in the superpoten- 
tial, see eq. @, leads to mixing of the neutralino and neutrino states ||. We 



treat the effects of R-parity violation perturbatively, using the method in |25 
in order to diagonalize the various mass matrices and hence determine the 
corresponding couplings. In doing this we follow the conventions given in ref. 
[ p6fl . In two component notation let denote the column vector of neutrinos 
and neutralinos in the basis 

< = (Vl^LVl, -iA', -iA 3 ,Vk^ 2 ), ( 51 ) 
where are the neutrino fields, —iX' and — i\ s are the unmixed photino and 

1 2 

gaugino states, respectively, and if)jj are the neutral components of the Higgs 
doublets. The Lagrangian describing the neutrino/neutralino masses is then 
given by, 

£ maS s = ~% T M % + h.c. (52) 
The mass matrix can be written in a block form, 
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m 



M ° = { m - m x , ) ■ < 53 > 

Here, the submatrix m contains entries from the bilinear Jft p parameters, 



m 



■\gfur 







hgw, 

{gu^ 
\gu T 



ijj; 



(ui). M. x o is the MSSM neutralino mass matrix, given by, 



M x o 



Mi 


-\g'vi 
V \g'v2 



o 

M 2 



\gvi 
-\gv2 



-\g'vx 

o 

-n 



yv 2 \ 



-2^2 





One defines the mass eigenstates by means of the rotation 



0,j 



The mass matrix then is diagonalize by 

E*M<0 = M , 



(54) 



(55) 



(56) 



(57) 



where A4q is a diagonal matrix, with the eigenvalues as elements. Note that 
these conventions are chosen such that in the limit of vanishing lfi p parameters 
the usual MSSM convention ilfl, H — »■ N ', is recovered. 



7.2 Charged fermion mass matrix 

For the charged fermion sector the mass part of the Lagrangian can be written 

as, 



EL* = -*£)-M±*( +) + h.c. 



where in two component Weyl spinor basis 

The (5 x 5) charged fermion mass can be written as, 



(58) 

(59) 
(60) 
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M 



± 



Mi E 
E' M x ± 



where M x ± is the usual MSSM chargino mass matrix 



M x ± 



(61) 



(62) 



v 75^ 2 » ) 

and the sub-matrices E and E' give rise to chargino-charged lepton mixing. 
They are defined as 



E 



E' 



( -j^9u e e e \ 



-j^hiiLUi -^h 2 iUJi ^h 3 iUJi 



(63) 



(64) 



Mi finally is the charged lepton mass matrix due to the ordinary Yukawa 
interactions. The rotation from the weak to the mass eigenstate basis is defined 

as, 



(±)i 



(65) 



where A^ diagonalize the mass matrix via 



(A-)*M±(A + y = Diag{m l i ,m+} 



(66) 



ml are the known masses of the charged leptons, and m ± the "chargino" 
masses, k — 1, 4. 

Note, that without making assumptions about the structure of Mj, it is 
impossible to diagonalize A4±. For numerical purposes one therefore has to 
assume that Mi (i.e. the Yukawa couplings h^) can be brought into a diagonal 
form. The constraint that the smallest 3 eigenvalues of Ai± should coincide 
with the known masses of the charged leptons, can than be reformulated as 
a set of equations for the Yukawa couplings, which for finite Ifl p parameters 
differ from the MSSM equations. 



16 



8 Appendix B: Approximate diagonalization 
of mass matrices 



8.1 Neutral fermion mass matrix 

If the 1/lp parameters are small in the sense that for 

£ = m ■ M-J (67) 

all <C 1, one can find an approximate solution for the neutrino/neutralino 
mass matrix H. Explicitly, the ^ are given by 

g'M 2f i 

61 = 2det{M x «) K m 
^ 2 = -2det(M x0 ) K m 

^ ~ m + Adet(M x o) A " (70) 

(ggjfr + g' 2 M 2 ) Vl 
Adet(M x o) 

It is interesting to note, that only £ i3 is not proportional to Aj, where Aj is 

Aj = iiuJi - vi€i. (72) 
In leading order in £ the mixing matrix S is given by, 



vj o \( i - |ee f -e 
o iv* e f i - |e + e 

#T iv*(i-|^) I 



(73) 



The second matrix in eq. ([T3|) above block-diagonalize Aio approximately to 
the form diag(m e //, A4 x o), where 



! T Miff 2 + M 2 fl' 2 



/ A 2 A e A M A e A r 

A e A M A 2 A M A T |. (74) 
V A e A T A M A T A 2 
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Here, det(A4 x o) is the determinant of -M x o 

det{M x o) = -^ 2 M X M 2 + -^v 1 v 2 (M 1 g 2 + M 2 g' 2 ). (75) 

The submatrices N and V v in eq. (173) diagonalize .M x o and m e ff in the 
following way: 

N*M x oN j = diag(m x o), (76) 

Vjm e ffV v = diag(0 } 0, m„), (77) 

where m v o are the "neutralino" mass eigenstates, and the only non-zero neu- 
trino mass is given by 

m„ = Tr(m eff )= ^ \A\ . (78) 



Note that this convention is chosen such that eq. (|76| ) coincides with the 
standard MSSM convention by Haber and Kane M. 
For V v one can find an analytic form, 



/ 1 \ / cos # i3 - sin 13 



i.e., 



cos# 23 -sin 6 23 ■ 1 | . (79) 
y sin #23 cos #23 / \ sin #13 cos 6*13 

(cos #13 — sin #13 \ 

- sin #i 3 sin #23 cos #23 - cos #13 sin # 23 , (80) 
sin #13 cos #23 sin # 23 COS#i 3 COS# 2 3 / 

where the mixing angles can be expressed in terms of the "alignment vector" 
A as follows: 

tan #13 = ' 1 1 



(A£ + A*)s 
A, 



tan #23 = -f-. (82) 
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8.2 Charged fermion mass matrix 

Similarly to the case of the neutralino mass matrix discussed above, for small 
values of the ]fi p breaking parameters it is possible to find an approximate 
diagonalization procedure for the charged mass matrix. 

Using the fact that the matrix E' is proportional to the Yukawa couplings 
of the charged leptons, one can assume E' w 0. 

Define then, 



a = e{m x ± 



C R = M}EM-UM~l) T = M}e L {M-i) T . 



Note that £r ~ £ L (mi/M). Explicitly. 



(f*). 2 = !i g 2 v2^i 

[i 2fi(Mfi — \g 2 viv 2 ) 
The rotation matrices A are then given in leading order in £^ and £r as 



o u* j\ e L i - uui 

(A+)t= ' 



The conventions are chosen such, that in the limit of vanishing I/lp eqs 
simplify to the MSSM notations: 

U*M X ±V ] = M x ±, 
where again the hat indicates a diagonal matrix. 



83) 
84) 

85) 
86) 

87) 



89) 
90) 
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V V 

Figure 1: Diagram illustrating the connection between the Majorana mass of 
the neutrino and the amplitude of double beta decay. 
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Figure 2: Decomposition of the double beta decay amplitude into long-range 
and short-range contributions. 
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Figure 3: Leading double-beta- decay Feynman graphs in the Ifl p model. 
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Figure 4. a: Upper limit on the individual parameter t\ for tan/3 
uji = as a function of (M 2) iJi). 
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Figure 4.b: Same as figure J^.a, but for oj\, assuming e\ = 0. 
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Figure 5: Calculated half-life for the (3[3q u decay of 76 Ge as function of uj\ for 
a random variation of the MSSM parameters, M2 and fi from 100 GeV to 1 
TeV and tan/5 =1 — 50. 
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Figure 6. a: Upper limits on (ei,ivi) assuming (m u ) < 0.5 eV, tan/? = 1 and 
H = 100 GeF ; for different values of M 2 , M 2 = 100,200,500, 1000 GeV. Note 
that the allowed range is always in between two lines of constant M 2 . 
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Figure 6.b: Same as figure 6. a, but for tan (5 = 10. 
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